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ABSTRACT  :  A  set  of  simultaneous  ordinary  differential 
equations  is  derived  for  the  peak  pressure  and  time 
factor  of  shockwaves.  These  equations  satisfy  the 
Ranklne-IIugonlot  conditions  ancf  the  partial  differential 
equations  of  fluid  dynamics.  No  particular  assumptions 
are  made  as  to  the  nature  of  the  medium,  so  the  equations 
hold  for  any  gaseous  or  liquid  medium.  For  application 
to  explosion  phenomena,  expressions  for  the  initial 
conditions  are  found  for  the  Interface  between  the  explo¬ 
sion  products  and  the  surrounding  medium  (which  in  this 
case  is  the  place  where  the  shockwave  is  formed).  Rela¬ 
tionships  axe  derived  for  tvo  cases : 

(1)  that  the  shockwave  is  observed  at  a  fixed  point 
(applicable  to  the  measurement  of  airblast 
waves) . 

(2)  that  the  point  of  observation  moves  with  the 
medium  (applicable  to  shockwave  measurements 
in  water) . 
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A  THEORY  OF  THE  PROPAGATION  OF  SHOCKWAVES 
AND  THEIR  FORMATION  BY  EXPLOSIONS 


I  INTRODUCTION 


1.  It  has  long  been  recognized  that  the  disturbances 
emitted  from  explosions  are  shockwaves,  l.e.,  waves 
having  an  infinitely  steep  front.  Such  shock- fronts  are 
fully  described  by  the  well  known  Ranklne-Hugoniot 
conditions,  whereas  for  the  phenomena  behind  the  front 
the  hydrodynamic  equations  of  an  inviscid  fluid,  neglecting 
gravitational  effects,  are  applicable.  It  is  the  objective 
of  this  paper  to  combine  these  relations  in  such  a  way 
that  information  on  the  change  of  the  pressure  of  the  front 
with  distance  can  be  obtained.  This  can  be  dons  only  if 
the  distribution  of  the  pressure  behind  the  front  is  taken 
into  account.  In  our  case,  the  pressure  distribution 
immediately  behind  the  front  is  sufficiently  described  by 
two  parameters,  the  time  factor  ci  and  the  shape  factor 
cft/c^.  In  this  paper  two  different  definitions  for  the 
time  factor  are  used,  namely: 


(a.) 

and 

(G) 


CK 


/■  V  **J  c. 


4  xc 


*. 


c. 


P  and  its  derivatives  refer  to  the  pressure  immediately 
behind  the  shock-front;  £  is  the  time,  f  the  radial  distance, 
^  the  velocity,  a,  a  reference  length  (in  the  study  of 
explosion  phenomena,  the  radius  of  the  unexploded  charge) 
and  a  reference  velocity  (for  instance,  the  sound 
velocity  of  the  undisturbed  medium).  These  two  reference 
magnitudes  are  introduced  in  order  to  obtain  the  time 
factor  as  a  dimensionless  magnitude.  The  time  factors 
o(  and  c((u'  are  related  to  the  initial  decay  of  the  pressure- 
time  curve  which,  in  the  case  of  c( 
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Figure  1 

SKetch  of  a  pressure- time 


diagram. 


is  obtained  by  observing  the  pressure  from  a  fixed  point, 
but  in  the  case  of  c(vu^ ,  from  a  point  moving  with  the 
medium.  These  two  cases  have  analogies  in  the  experimental 
measurement  of  shockwaves  due  to  explosions.  In  air  blast 
measurements,  the  pressure  recording  gauges  are  usually 
so  rigidly  fixed  that  they  cannot  move  significantly  under 
the  action  of  the  impinging  blast  wave.  Here  ot  applies. 

For  the  measurement  of  underwater  shockwaves,  tne  gauges 
are  usually  not  very  rigidly  mounted.  Considering  the 
high  pressures  of  the  shockwave  and  the  relatively  small 
particle  displacements  in  water,  the  case  o('u'  seems  to  be 
more  appropriate  for  underwater  explosions.  However,  it 
should  be  noted  that  only  small  differences  are  to  be 
expected  from  these  two  treatments  in  the  range  practical 
shockwave  measurements. 


2.  The  analysis  presented  in  this  paper  yields  two  simul¬ 
taneous  ordinary  differential  equations  for  the  peak  pressure 
and  the  time  factor  as  functions  of  distance.  In  the 
equations  for  the  latter  magnitude  the  shape  factor  appears 

as 


or 


a 

of* 


These  magnitudes  are  related  to  the  second  derivative  of 
the  pressure  with  respect  to  time  immediately  behind  the 
front.  Again,  the  first  of  these  refers  to  a  fixed  point 
of  observation,  the  second,  to  a  point  moving  with  the 
medium. 
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3.  To  integrate  the  differential  equations,  the  shape 
factor  must  be  known  as  a  function  of  distance.  This 
magnitude  depends  only  on  the  shape  of  the  wave;  no 
information  about  tho  pressure  or  the  time  scale  is 
necessary  for  its  evaluation.  For  Instance,  if  the 
pressure-time  curve  (observed  from  cither  of  the  points 
of  reference  discussed  above)  is  an  exponential  curve, 
the  shape  factor  has  the  value  one;  if  it  is  a  straight 
line  the  shape  factor  is  zero.  Shockwaves  from  explosive 
sources  have  the  important  property  that  their  shape 
factors  are  constant  within  an  appreciable  range  of  dis¬ 
tance.  This  holds  true  particularly  for  the  initial  portion 
of  an  underwater  shockwave  which,  at  all  distances  experi¬ 
mentally  investigated,  is  an  exponential  curve  to  a  surpris¬ 
ing  degree  of  accuracy. 

This  fact  has  first  been  used  by  Kirkwood  and  Bethe  [a]* 
in  the  calculation  of  shockwave  phenomena.  Following  their 
usage,  we  call  this  procedure  the  "peak  approximation". 

In  our  analysis  -  which  is  rigorous  up  to  this  point  -  this 
amounts  to  an  a  priori  assignment  of  a  value  to  the  shape 
factor. 

4.  No  assumption  of  a  particular  shape  of  the  wave  is 
incorporated  in  the  equations  derived  here.  This,  together 
with  the  rigorous  character  of  these  equations,  makes  it 
possible  to  study  the  validity  of  the  peak  approximation 
and  to  show  hovr  sensitive  the  solutions  are  to  the  assumed 
shape  of  the  wave. 

5.  In  order  to  obtain  numerical  solutions  of  the  differ¬ 
ential  equations,  the  initial  conditions  must  be  known. 

In  the  calculation  of  shockwaves  from  explosions  these 
conditions  are  provided  by  studying  the  phenomena  at  the 
interface  between  the  reaction  products  and  the  surrounding 
medium,  l.e.,  by  studying  the  formation  of  the  shockwave. 

It  has  often  been  pointed  out  that  a  high- amplitude  pressure 
pulse  must,  during  propagation,  change  its  shape  in  such  a 
way  that  a  shock-front  is  finally  built  up.  This  phenomenon 
seldom  occurs  with  shockwaves  due  to  explosions.  The 
rapidity  with  which  the  pressure  is  built  up  by  an  explosion 
causes  the  pressure  wave  to  have  a  steep  front  from  the  very 
beginning.  It  is  a  shock- front  which  emerges  from  the 
surface  of  explosion  products  into  the  surrounding  medium 


*  All  such  letters  refer  to  the  list  of  references  at 
the  end  of  this  report. 
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a  very  short  time  after  the  detonation  of  the  charge. 
Therefore,  this  front  can  provide  the  initial  conditions 
for  the  differential  equations  for  the  shockwave  peak 
pressure  and  for  the  time  factor. 

6.  Several  similar  approaches  to  the  problem  of  shock¬ 
waves  due  to  explosions  have  been  made  previously,  notably 
by  Kirkwood  and  co-workers.  The  Kirkwood- Be the  theory  [a  ]  - 
though  very  successful  in  its  application  to  underwater 
explosions  -  is  based  on  a  propagation  theory  which  is  not. 

a  rigorous  solution  of  the  hydrodynamic  equations.  The 
Kirkwood  Brinkley  theory  [b  ]  which  is  applicable  to  any 
fluid  is  a  rigorous  approach  based  on  energy  considerations 
and  a  hybrid  form  of  the  Euler  and  Lagrange  equations  of 
fluid  dynamics.  The  peak  approximation  (or  "similarity 
restraint"  as  it  is  called  in  that  paper)  is  incorporated 
in  the  final  equations  by  assuming  an  exponential  wave. 

7.  In  the  present  paper,  an  attempt  is  made  to  treat  this 
problem  in  a  straightforward  manner,  starting  with  the 
Euler  equations.  A  similar  approach  has  been  made  previ¬ 
ously  in  other  papers  [ c  —  f ]  In  these  papers  equation 
1.6  or  similar  expressions  are  derived  in  various  ways, 
but  the  treatment  is  not  carried  much  further. 


II  HYDRODYNAMIC  RELATIONS 

8.  The  fundamental  equations  for  the  fluid  motion  of  an 
inviscid  nodium,  neglecting  the  influence  of  gravity,  are 


(°l) 

^  + 

-CC  ACp, 

+  -  0 

2-u/> 

+  ■**  —f~  ~  0 

(02) 

A  + 

^/v 

M 

4JL 

-  0 

The  subscripts t  or  r denote  the  partial  derivatives  v 1th 
respect  to  these  magnitudes.  In  general, 


Throughout  this  paper, P  denotes  the  pressure  in  excess  of 
the  pressure  of  the  undisturbed  medium,  . 
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In  the  above  equations  ^  is  the  particle  velocity,/? 
the  density  and  S  the  entropy.  and  £  are  the  space  and 
time  coordinates  respectively.  Equation  (0.3)  states  that 
we  are  considering  reversible  processes,  which  Is  consistent 
with  the  assumption  of  an  inviscid  fluid.  Equatic.iS  (0.1) 
to  (0.3)  are  not  applicable  across  discontinuities  1. shock- 
fronts;,  as  discussed  below. 

9.  Equation  (0.2)  is  written  for  the  case  of  spherical 
symmetry.  The  case  of • cylindrical  symmetry  is  obtained 
by  omitting  the  factor  2  of  the  last  term  of  (0.2),  the 
plane  wave  case,  by  omitting  the  last  term  entirely. 

10.  Since 


-  /(/  +  -A  i  $) 

the  density  can  be  eliminated  in  (0.2).  Introducing  the 
sound  velocity 


9  -*->.)  ] 

V 


We  obtain  from  (0.2)  and  (0.3) 


(oh) 


A  * 


+ 


D 


11.  Introducing  the  notation 


i dt  ^  Dr 


(0.1)  and  (0.4)  can  be  combined  to  take  the  form 


(os) 

(O') 

/  (c-m) 

- Jt>  —  X* , 

A 

2c U 

r 


2cU 


•y^’^can  be  interpreted  as  the  derivative  of  with 
respect  to  time  along  a  "(c ,  +  m.)  -  characteristic",  i.e. 
along  a  curve  for  which  in  the  r-t  plane 


jL  t 

-  -  G  +  IV 

d-t 
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(f’) 

In  general,  ft  denotes  the  derivative  along  a  g- charac¬ 
teristic  for  which 


12.  If  a  fully  isentropi*  fxow  is  considered  (cSt  -  0,  ^  =  0, 
therefore  also^e^u^=  0).  we  may  introduce  the  Riemann 
function  which  is  defined  by 


(0.5)  and  (0.6)  can  then  be  transformed  into 


m 

(or) 


/  ^  \(t 

2ajlc 

(<r+  - 

r 

2ajo 

\° 

r 

which  is  the  well-known  Riemann  formulation  of  high- amplitude 
wa^es.  These  equations  are  not  always  applicable  in  the 
cases  considered  here. 


13.  For  the  shock- front,  these  equations  are  supplemented 
by  the  Rankine-Hugoniot  conditions.  These  are  three 
equations  which  permit  .<*,/>  ,  0  and  the  propagation  velocity 
of  the  front. u  ,  to  bo  expressed  in  terms  of  ^  ,  as  soon 
as  the  equation  of  state  of  the  medium  considered  is  known. 
Here ,  denotes  the  amplitude  of  the  shock- front.  There 
are  tables  available  which  give  these  relations  for  air, 
sea  water  and  fresh  water  [g  -  k].  For  our  purpose  we 
can  consider  that^,  p  ,  c,  and  U  are  known  functions  of 
at  the  shock-front. 


Since  we  are  interested  in  the  change  of  the  shockwave 
peak-pressure  (i.e.,  the  front  pressure)  with  distance, 
we  use  the  notation 


dJ 


•D  f  T>f 

—  +  u  — - 

'it  'SX 
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That  is,  we  arbitrarily  define  the  derivative  along  the 
u  -  characteristic  (which  is  equivalent  to  the  shock- 
front)  as  the  total  derivative. 

l4„  The  concept  of  the  characteristic  gives  us  a  geo¬ 
metrical  approach  to  our  problem.  Consider  the  ‘character¬ 
istic  triangle"  below. 


The  lines  in  this  figure  are  drawn  in  at-i*  plane  in  such 
a  way  that 


along 

AB 

AE 

AF 

FB 


d,t 
aL> “ 

d-t 

air 

<xt 
at  r 

oit 


/ 

c 

/ 

My 

I 

u~ 


We  then  call 


AB 

the 

•  I  G»  +  \JLf  — 

characteristic" 

AE 

the 

"  C.-LO  - 

characteristic" 

AF 

the 

"  tu  - 

characteristic" 

FIEB 

the 

"  U  - 

characteristic" 
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These  lines  are  not  straight  lines,  since,  e.g.,  1/U  at 
E  differs  generally  from  the  1/u  at  B.  This  character¬ 
istic  triangle  is  based  on  the  fact  that  U  which 

holds  generally  and  accounts  for  shock  phenomena. 

15.  In  Figure  2,  the  line  FIEB  represents  the  shock- front. 

The  regime  below  this  line  is  that  of  the  undisturbed 
medium  into  which  the  shock  propagates;  the  regime  above 
represents  the  shockwave,  where  there  are  high  pressures 
and  a  particle  velocity  towards  the  front. 

16.  Considering  the  equations  above,  we  find  that  (0.5) 
holds  along  the  line  between  A  and  B  and  (0.6),  between 
A  and  E.  Furthermore,  along  FIEB  the  Rankine-Hugoniot 
conditions  hold  and  along  FA  equation  (0.3)  holds.  Thus, 
at  A  we  have  the  same  entropy  as  at  F,  whereas  at  B,  in 
general,  the  entropy  is  different.  This  is  because  shock- 
fronts,  which  involve  irreversible  processes,  cause  an 
entropy  increment  which  depends  on  the  amplitude  of  the 
front. 

17.  As  far  as  the  thermodynamic  state  is  concerned,  we 

note  that  along  FA  the  common  adiabatic  (i.e.,  the  isentropic) 
holds,  whereas  along  FB  the  Hugoniot  adiabatic  applies. 

Since  both  of  these  thermodynamic  relationships  are  known, 
we  can  express  the  thermodynamic  state  along  FA  as  well  as 
along  FB  in  terms  of  the  pressure  along  these  lines,  if 
the  state  at  F  is  known.  The  same  holds  for  the  particle 
velocity^  along  FB,  but  not  along  FA. 


Case  I  Fixed  Point  of  Observation 

18.  The  Peak  Pressure.  Using  the  characteristic  triangle, 
we  will  now  derive  the  differential  equation  for  the  shock¬ 
wave  peak  pressure  as  a  function  of  distance.  In  Case  I, 
we  consider  the  point  of  observation  as  fixed. 

On  integrating  (0.5)  and  (0.6)  along  AB  and  AE  respec¬ 
tively,  we  have 


CO 

02) 


(ft ij-**  +  u< 

£77£(VA)-(*A-“f) 
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vhere  the  double  subscript  denotes  the  average  obtained 
In  the  Integration. 


In  this  paper,  we  will  be  primarily  concerned  with 
first  or  second  order  expressions.  Hence,  these  averages 
need  to  be  correct  only  to  the  first  order.  For  Instance, 
//sc  then  becomes: 


db  *  •  I  1  /  \ 


With  this  approximation,  (1.1)  and  (1.2)  can  be  combined 
to  give 


((3) 


■  -  [(vi  +  Pr)]  <*-*)■ *[(7). +  (t\]  (n) 

19.  If  we  expand  (1,3)  by  means  of  a  Taylor's  series, 
terms  of  order  zero,  one,  two,  and  higher  appear.  Consid¬ 
ering  only  first  order  terms,  equation  (1.3)  becomes 


— —  (*» -O  -  —  *<(**-%)-- 

p  C  dbt  '  •  />c  '  t  v  *  1 ' 


/ 


/c 


+  -f-  U- 1 )  +■  o  — /t  -  t.)-o 

aif,  ait  l  «  */  >-  V  »  V  h-  V  *  £/ 


All  magnitudes  without  subscript  refer  to  the  point  I, 
i.e.,  to  a  and  £  .  For  the  various  differences  in^  , 
we  obtain  the  following  first  order  relations: 


t  -  t 

A  1 


C  r_ 

U 

t  +  Ur 
C  4-  Lt  17 

u)‘ 

2  c  £  r 

(t  ♦ 

gr 

(<■  *■  cl)  (c *U-u) 
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Substituting  these  expressions  in  (1.4)  and  dividing  by 
4  r  ,  we  obtain  on  collecting  terms: 


('■'>  -u  % 


Using  the  notation 


and 

y*  _ 

(1.5)  becomes  ° 


F«a'° 

where  P„  and  Pa  are  functions  of  the  pressure  alone : 


When  a  double  subscript  appears  in  functions  such  as  these, 
the  first  subscript  refers  to  the  case  being  discussed 
(I  or  II)  and  the  second  subscript  to  the  number  of  the 
function.  Thus  is  the  second  P- function  of  Case  I. 

20.  The  differential  equation  Just  formed  has  been  derived 
by  other  methods  (see  paragraph  7)*  For  instance,  equations 
(0.1)  and  (0.2)  have  been  combined  to  give  (1.6)  directly, 
[d,  e]  . 

21.  Th<*  Time  Factor.  In  order  to  obtain  a  differential 
equation  involving  Iq/et*  ,  we  have  to  consider  the  second 
order  terms  in  the  expansion  of  (1.3) •  Rather  than  follow¬ 
ing  the  same  lines  as  in  the  foregoing  paragraphs,  we  will 
obtain  the  desired  equation  by  an  analytic  method  in  which 
the  computations  are  shorter.  The  two  procedures  are 
analogous,  however. 
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Differentiating  (0.2)  and  (0.4)  partially  with  respect 
to  ,  ve  have 

(2/)  +•  u.  u4  a,.  y  A  ^  0 

(2  2)  At  +  ^  ^  +  %  4.  4  />**  +  “v  0‘Vt  +  “T  (f  ''I 

4.  i£L 

Differentiating  (1.6)  along  U  yields 

»  A  A,  iS.)  *  «./„  .  P  45.) 

^  {  dp/  «  l  "  ■O'  "  dp  S 


-  0 


(«> 


<*fi 

i  „  dp„  _ 

*  a  7*  '  77 


In  addition, 


d&  y6  7 d^>  ^ 


ctt* 


dp  dii  dp?  \d* 

& 

so  that  we  obtain 

«rf£\  . .*  <*"• 


2  da.  d-51 

-  u  — v  +  O 

dj) 


<dH 


[_dp  dp 


U 


<*P*J 


) 


u 


f  dUi  di)  d*u.  ^ 

■"  ^  ^  * u  ^7 


df  vr  0^  ”  */>, 

eta,  d 
pLfc  djfi 


ft7  r  7  du,  dVi  »  f  cLu  d(J 

a[  [_u  P*  ~7f  Jj,  +  u  P"  wf  TjT  +  J  u+V_ 


jj' cht_  cl(X 
df>  dr 


Moreover, 

(2S-)  d  /d“ 


—  (——)  -  Uu  -f  Jj-  /(J«  ) 

<u  \  ctl  J  **  * K  dt  '  *•/ 

Introducing  o(  as  well  as  (1.6)  into  (0.4),  we  obtain 
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where  q,0l  and  Q/ot  are  functions  of  -A'  alone.  Many  such 
functions  will  appear.  They  are  listed  for  ready  reference 
at  the  end  of  this  section.  As  above,  each  function  is 
identified  by  a  double  subscript:  the  first  number,  1, 
refers  to  Case  I,  and  the  other  two  give  the  sequence 
number  of  the  function. 


22.  By  differentiation, 


( u ?  (UQ'«  *  -*%r, 

( uPi • 


Combining  this  equation  with  (2.4)  and  (2.?),  we  obtain  an 
expression  for  as  follows: 


k  -  2isl  <**-%-  (X  + 


-  Utt, 


Substituting  for in  (2.1).  we  may  eliminate  between 
this  equation  and  (2.2)  to  obtain  a  single  equation  in  which 
Att  and  -Ptr  appear.  We  eliminate  by  using  the  follow- 
ing  relation: 

+»  +  u+‘*  ~  ~ 


so  that  we  obtain  the  equation 
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23.  Solving  for  ju.  ,  ^  and  p  ,  we  obtain  the 

following  results:  t  rr  t  t 


it  — 

<?/0J  0 

\A  - 

~hT  + 

u  ** 

9m  -t- 

& 

* 

r 

^  tr 

-  5.  ♦ 

/**. 

TV 

\ 

2 

C.  “ 

-  £/** 

Sjls.  Q( 

* 

K 

*. 

A  = 

1 

1 

iuz.  a 

a 


Further,  we  define 
Substituting  in  (2.6) 


yields  the  following  equation: 


(Z7) 

or 


era 

dr 


+ 


0 


Id- 


+ 


a** 


dOt 

otX 


rn 


3l« 


+ 


^5  & 


^/-i-  ^r)  -  0 


Again,  the  functions  denoted  by  P  depend  on  the  shockwave 
peak  pressure  only. 

2*+.  The  functions  denoted  above  by  P  and  Q  terms  are 
explicitly: 


Q101 

Q/ot 

Qic.i 


+  c^ju-u) 

1  t 


u 


* 


UP, 


—  ♦  u 


f  U0«+  1^, 
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Case  II  Moving  Point  of  Observation 

25.  The  Peak  Pressure.  In  Case  II,  it  is  assumed  that  the 
point  of  observation  moves  with  the  particles  of  the  fluid. 
The  approach  to  a  differential  equation  for  the  shockwave 
peak  pressure  is  much  the  same  as  the  one  given  in  para¬ 
graph  19  for  equation  (1.6).  The  characteristic  triangle, 
Figure  2,  takes  a  slightly  different  form,  as  the  line  AF 
must  now  be  considered  Instead  of  the  line  AI.  The  initial 
steps  are  the  same,  so  that  we  arrive  at  equation  (1.3) 
again.  Introducing 

*  A  *  »A 


instead  of ,  we  obtain  on  expanding  (1.3) 


-h.  (t-  t\  2  M  & 

fit  m  r»  cr)  ~  r 
(3  /)  '  r 

oLj,  cLt  V  «  »/  ¥  \ 

%  fa 

_  t  'N  ?Cu  (i  t  \  -  n 

e  W  y  *r)  ~  0 

where 

te~tr  = 

u.  ¥  t  *  U  (  A  ^ 

±  A.  _ 

/t  _t\ 

\  tf  ~ 

t  ♦  U  -  U.  (  A  v 

2a  (U  -u.)  /  \ 

f8 

(au-L))(tf  U-u)  \  A  v 

/t  M  \ 

u.  ♦  t-  U  (  * 

Making  these  substitutions  in 
by  ±A  -t f  we  obtain 

(3.1)  and  dividing  through 
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For  consistency  of  notation,  let 

/vH  ^ 


so  that  (3.2)  becomes 

vhere 


+  p„  =  0 


o  -  V  **  c3  (U  ~  u) 

21  "  (V-aH- 

p  -  _ >ea  IV"  (u-  u.)?1 

??  u  c*  + 

26.  The  Time  Factor.  The  procedure  of  obtaining  a  differ¬ 
ential  equation  for  the  time  factor  parallels  that  for 
Case  I. 

Here  we  define  qf^  by 

<  -  f&i&T 

<3*)  \ 

=*  (+**  ■*  +  +  “hr  +  '“V  +r) 

The  differential  equation  has  the  form 

(31)  —■  +  ^  +  F»  «N*  ('  +  /■»  -$br)  ”  ° 

where  the  ^  's  can  be  expressed  in  terms  of  the  functions 
and  ^  in  an  analogous  way  to  Case  I: 


1  +-  pn 
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4>tr  " 
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III  THE  FORMATION  OF  SHOCKWAVES  BY  EXPLOSIONS 

27.  The  Initial  conditions  for  the  differential  equations 
derived  above  can  be  obtained  for  the  case  of  explosion 
phenomena  by  considering  the  motion  of  the  interface  be¬ 
tween  the  gaseous  reaction  products  and  the  surrounding 
medium.  The  conditions  of  continuity  in  pressure  and 
particle  velocity  yield  all  the  necessary  equations  for 
the  initial  pressure.  The  use  of  the  characteristic 
triangle  is  particularly  convenient  for  determining  the 
initial  time  factor. 

Initial  Conditions  for  the  Steady  State  Explosion 

28.  Let  us  consider  the  idealized  steady  state  explosion.* 

It  is  assumed  that  the  explosive  is  converted  to  gaseous 
reaction  products  in  such  a  way  that  at  the  very  first 
moment  the  pressure  is  constant  throughout  the  space 
previously  occupied  by  the  explosive,  and  all  particles 
are  at  rest.  For  the  very  beginning  of  the  expansion  of 
the  reaction  products  at  the  boundary  we  have  the  conditions 

&0  u.  -a,  4’,  '  -&i 

(st)  '  -  07* 

The  asterisk  refers  to  the  reaction  products  of  the 
explosive,  the  symbols  without  asterisk,  to  the  surrounding 
medium.  Equation  (5*1)  represents  the  condition  of  con¬ 
tinuity  of  pressure  and  velocity  at  the  interface.  Equation 
(5.2)  determines  the  velocity  resulting  from  the  expansion 
of  the  initially  motionless  reaction  products.  Since  these 
(at  least  in  the  moments  we  are  considering  here)  undergo 
only  expansions,  which  are  necessarily  isentropic  -  the 
Riemann  function  is  used: 


h 


♦It  would  be  more  realistic  to  assume  that  the  charge  is 
initiated  at  the  center  and  that  a  spherical  detonation 
wave  spreads  through  the  charge  causing  the  explosion  of 
the  particles  in  concentric  shells,  one  after  another. 
This  case  can  be  treated  approximately.  It  will  be  pre¬ 
sented  in  a  later  publication. 
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Equations (5.1)  and  (5.2),  together  with  the  equations  of 
state  and  the  Rankine-Hugoniot  equations  determine  the 
pressure  and  the  velocity  .  The  procedure  is  well 
known  and  will  not  be  repeated  here :  see  for  instance  [ 1 ] - 
This  yields  the  initial  conditions  for  the  differential 
equations,  giving  the  pressure  versus  distance,  (1.6)or  (2.3). 

29.  Case  1.  We  will  now  derive  an  expression  for  the 
initial  time  factor. 


Characteristic  Triangles  at  the  Interface  between 
Reaction  Products  and  Medium  for  the  Steady 
State  Explosion 

In  Figure  3»  the  line  OelK  represents  the  interface  be¬ 
tween  the  reaction  products  (left  hand  side)  and  the 
surrounding  medium  (right  hand  side).  The  regime  below 
feO  is  that  of  the  unreacted  explosive.  Between  aef  we 
have  the  reaction  products  at  steady  state  conditions. 
The  expansion  takes  place  in  the  regime  above  Klea.  The 
line  1EB  represents  the  front  of  the  shock  which  spreads 
into  the  surrounding  medium.  At  el  the  pressuie  of  the 
reaction  products  drops  instantaneously  f rom  ^  to  ^  . 

30.  To  summarize,  we  have 


(S3)  - 0  0  “l  0  ° 


(Si)  u,“  t  0  -  u*  «•  <£  --  ul  *  az 


21 


NAVORD  Report  2195 


Following  the  same  lines  as  above,  we  obtain 
(JS)  -(«.  -*«)  *  ('•-*•) 

(*i  ■  «i) +  -  O  • 

With  (5A)  and  ,  we  find  by  adding  the 

quantities  in  (5.5) 

^)„  (*<  ~  *')  *  &  ~  *?*)  +  "r- 
Considering  the  first  order  terms  only,  we  obtain 

^  ft  */?)  ^  ‘.)  ♦  (£  - fr- 

(v,)-tVv^ 


?u.C 


All  magnitudes  without  subscripts  refer  to  the  thermo¬ 
dynamic  state  at  the  point  1. 

30.  The  differences  in  t  can  b?  expressed  as  follows: 


/  \  U  -tu 

*«  -  *•  "  (*«  _t') 

*«-%  -  !-(*<-*■) 

__  C.*- 

*b  ta-  ”  t|f  _t'  3  c  </*-•- u-*+ c-* 
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so  that  (5*6)  becomes 


(s?)  (—  4 

V/e  /’t/  * 


.  _  _L )  / _ s^_\  +  1± t 

\  dp  jfcj(c.+  U-oL//  oit  r 


can  be  eliminated  by  means  of 
,  («)  u-h 


\  +  ia±  ( u~u-  \ 

/  dut  r  It  +  U-n-/ 

_  *  *r”  *  «  i  ~ 

2u.j  /  -  i  -  Cc  n 

j  *2u<\.Jd7) jj 


60  -  +  - 

dr-  U 


Further  substitution  for  cCjo/af^  and  ,^6t  yields 

/jT/J  -  (-L  +  lit  or,  -  P-["L/—  +  4-  e(^  - ' 

(•"/  l^e  /V/l  U  /  <le  '  r#[_l/>e  /t/ 


-  [a  A-  +  Vt)+  — «-  —P=^~) 

Ol  o  [_  yd  /  t  /  c.*U-^\o^b  /><VJ  ^  \c4U*a/ 

*  (  .8,«  -  a 

r0  \ 3  t e  +  c  / 

where  <T,  refers  to  the  initial  state  of  the  surrounding 
medium  at  ^  =  <2^  .  Solving  (5.8)  for  ^  ,  we  obtain 


or,  =  a, 


tit+Gn*??  Gn 

•'  6  i.^v 

L  *  /c  /  c  J 


where 


Q  *  — 
ii 


Gu  “ 


■P^O  (U  -  “•)  +  ^ 

2  u-c.(U-  a)  P„  ix  c.  U  f?i  /  <6^  ^ 

t+LJ-UL-  /»C-  OlJ-U-  Vo^  p  + 

2aV(2  c,VQ 
jC£  +  u.'  +-  c* 


_  /  W  U _ ^ 

t+U-U.  \  ot^>  /  l  /  \  f  ct  (u  -  u)  +  ul/  P/t  / 

G-functions  without  an  asterisk  refer  to  the  surrounding 
medium  and  depend  only  on  the  pressure^,  .  G{2  refers  to 

the  reaction  products  of  the  explosion. 
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31.  It  is  of  interest  to  note  that  if  we  had  set  up 
equations  (5*5)  using  the  characteristic  KB  instead  of 
KE,  the  result  would  have  been  the  same. 

32.  Case  II.  Retaining  in  (5.7),  we  obtain 

(b.O  _  Is  ,  +  A  -  -L)  4- i-)"|  -Mi  pa.  -LV_±_\ 

i  Zitif  U~n  \  2>gT c*  /  u  -  c  *  -  Cc*  71  n 

►.  (t*uu/  n  [_2  rT3t/7^%T*yJ  "  0 

Solving  for  ,  yields 


(6.2) 


2Ht(U-u.)  a  /Un  >  ^/  fc  \  .  ,  •  •  /  2t«*+ •«.*  \ 

(*)  /  ;-^rr-T. —  ~UPf/(~ar  ~  jr)K<-  +  u-Z)  +  *•  v3te.,V^«-~Iy 


+  '  )  AL. 

V*-  t  «■ '  “» 


^  U  /  du. _ M  f  c  \ 

a..  {  oLb  fc.J  (t+U-u-/ 


/v  11  _C«;  ^  f  G#/  +  GXi 

(6-V  a<  =  G«JJ^r 


where 


.  ( KJ- a) 


1 * U-  u- 


-  UR, 


--y-— ) 

/c  A  (J-u/ 


?  U-  C*  f  2  C  e  +  U*  *) 
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IV  PLANE  AND  CYLINDRICAL  WAVES 


33.  It  is  of  interest  to  know  how  the  functions  of  the 
foregoing  sections  would  look  in  the  case  of  plane  waves 
or  for  cylindrical  waves.  Once  the  spherical  case  has 
been  worked  out.  it  requires  only  very  simple  transfor¬ 
mations  to  obtain  the  equations  for  the  two  cases  mention¬ 
ed  above.  All  of  the  calculations  start  with  equations 
(0.5)  and  (0.6).  In  the  plane  wave  case,  the  term 
is  zero;  as  a  result  no  term  will  appear.  For  the 
cylindrical  wave,  the  divergence  term  is  -c-u^,  • 


Cylindrical 

(Ua.) 

(23  a.) 

(374 

(y./a) 

(s.fa) 

(*■30.) 


wave: 
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& 

4 

*  -  * 
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dp 
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itt°L  . o 

dnc 

2r 

O'. 

dOC 

t 

4 

/  +  «l)-o 

Vi*1 

2'  a/  V 

or 7 

dOf 


(*) 


4* 


Pjj_  +  '„<***.  P'S  a0*1 


a  -  -- 

7 


a*-  % 


'C  +  g*  *  —  a-1 
rt_  —  /r  111 


1  +  +  6 


-/»C  /  t 
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« -  G.\ 


—  ♦  t4  -6, 

v*  />*  ' 


(a) 


l  +  a** 


rr 


Plane  wave: 


Q  //) 

+ 

<Ly 

(23/) 

+ 

<*.r 

(3.7/) 

da  + 

(*/*) 

* 

<?,  <**('  +  Pnft)  =  0 

Pts  w001  (i  t  rtl  =  0 
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It  is  interesting  to  note  that  the  functions  Pt<<.  and 
are  the  sane  for  the  spherical,  cylindrical  and  plane 
wave.  Hence,  if  they  are  computed  for  a  particular  medium 
they  may  be  applied  to  any  one  of  those  types  of  waves. 

In  a  plane  wave  0(t  ,  Is  zero.  That  means  a  plane 
steady  state  explosion  produces  a  step  wave,  whose  pressure 
Is  constant  with  distance.  This  holds  only  fairly  close 
to  the  charge.  Later  the  rarefaction  of  the  explosion 
products  becomes  effective.  A  rarefaction  wave  follows 
this  step  wave  and  finally  overtakes  the  front  thus  causing 
the  pressure  to  drop.  These  phenomena  are  veil  known 
from  the  theory  of  the  shock  tube. 
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R.  H.  Matthias,  249  Kenville  Road, 

Buffalo  15,  N.Y. 

Prof.  P.  C.  Cross,  Dept,  of  Chemistry  and 
Chemical  Engineering.  University  of 
Washington,  Seattle  5,  Washington 
Prof.  Werner  Goldsmith,  Dept,  of  Engineering 
Design,  University  of  California, 

Berkeley  4,  o&iif. 
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